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1. INTRODUCTION. The method of completing squares yields an elemen- 
tary procedure to evaluate 
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Observe that Aac — b 2 > is required for the convergence of (fl]) . 

The goal of this paper is to present a new proof of ([2]). We illustrate a 
technique that will apply to any rational integrand. Providing new proofs of an 
elementary result, such as (J2J), is usually an effective tool to introduce students 
to more interesting Mathematics. The method discussed here has a rich history 
that we describe in section 2. 

It is an unfortunate fact that, despite our best efforts, evaluating definite 
integrals is not very much in fashion today. Thus we rephrase the previous 
evaluation as a question in dynamical systems: replace the parameters a, b, and 
c in ([!]) with new ones given by the rules 
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and to prove that 
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Once this is done, we can pass to the limit in Q and use the invariance of / to 
obtain 



7T 27T 
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lim a n ^Aac - b 2 ' 
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This leads directly to a proof of ©. The advantage of this method is that it 
generalizes to integrands of higher degree. 

We call (0) a rational Landen transformation. In section 2 we discuss the 
historical precedent and motivation behind such transformations. This history 
connects © to the magic of the arithmetic-geometric mean, 7r, and a wonderful 
numerical calculation of Gauss. 

The rest of the paper is devoted to a detailed proof of the Landen trans- 
formation: the invariance of the rational integral ((4]) and the evaluation of the 
limits in ([5]). A scaling of the integrand that is a crucial step in producing 
this transformation is presented in section 3. The following section presents the 
trigonometrical aspects of this problem and completes the proof of (HJ . An alge- 
braic calculation shows that the discriminant of the quadratic in ([1]) is preserved; 
that is, 

Aac-b 2 = iaia-bl. (7) 

This invariance is used in section 5 to analyze the dynamics of ([3]) and to 
establish ([5]). 

2. LANDEN TRANSFORMATIONS. Many of the evaluations encoun- 
tered in integral calculus illustrate the fact that definite integrals correspond 
to special values of functions. For example, the last integral in © is given by 
7r = tan _1 (oo) — tan _1 (— oo). Other special values appear in elementary courses: 

sin- (*) . (8) 



o V3 — x 2 \ \fi 

The same is true for more complicated integrals. For instance, when < 
b < a < 1, 

G{a,b):= = -K(k), (9) 

Jo V a 2 cos 2 9 + b 2 sin 2 6 a 

with k 2 = 1 — b 2 j a 2 . Here K is the complete elliptic integral of the first kind 
defined by 

K(k) = I' - dX (10) 

Jo ^{l-x 2 )(l-k 2 x 2 ) 

Elliptic integrals appear at the center of classical analysis. Their name comes 
from the fact that they provide explicit formulas for the length of an ellipse. 
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The inverse of 

dx 



m = 



o y/(l-x 2 )(l-k 2 x 2 ) 



is similar to sin z, so (j8]) and IT0|) are not so different after all. This new function 
is the elliptic sine (or sinus amplitudinus) of Jacobi [15], denoted by sn z. It 
completes the trilogy: sinz (circular), sinh z (hyperbolic), and sn z (elliptic). 
The question of evaluating definite integrals sometimes comes down to how 
many functions one knows. 

Our complaint that students today are exposed only to the most basic of 
functions is not new. Klein states [16, p. 294]: When I was a student, abelian 
functions were, as an effect of the Jacobian tradition, considered the uncontested 
summit of mathematics and each of us was ambitious to make progress in this 
field. And now? The younger generation hardly knows abelian functions^ 

Suppose that a and b are positive real numbers. It is not hard to check that 
the sequences {a n } and {b n } defined recursively by 

a-n+i = — ~) ^"+1 = yj a n b n , (11) 

ao = a, and bo = b converge to a common limit: namely, the arithmetic- 
geometric mean of a and b, denoted by AGM(a, b) . This is a fascinating function; 
the book [5] explains its connections with modern algorithms for the evaluation 
of 7r. The reader will find in [1] a survey of maps similar to (llip and an extensive 
bibliography. 

At the turn of the eighteenth century, Gauss [13] was interested in lemnis- 
cates and their lengths. After a numerical calculation, he observed that 

1 



AGM(1,V2) 

and 

2 f 1 dx 



7T Jo y/l-X A 

agree to eleven decimal places. (The integral gives the length of a lemniscate.) 
With remarkable insight, he discovered that the elliptic integral G(a, b) in ([9]) 
remains invariant if the parameters (a, b) are replaced with their arithmetic and 
geometric means; that is, 



G{a,b) = G\-^-,Vab\. (12) 

Iterating, passing to the limit, and using the invariance of the elliptic integral 
G yields 

= 2AG^ (13) 



The authors learned of this quote from the preface of [5] . 
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The convergence of the arithmetic-geometric mean iteration (fTTj) is quadratic, 
meaning that \a n+1 - AGM(a, b)\ < C|a„ - AGM(a, 6)| 2 for some C > 0. Thus 
(fTTj) leads to a rapid evaluation of the elliptic integral G(a,b). This iteration 
has been used for the numerical evaluation of elliptic integrals. See [7], [8], [9], 
[10], or [11] for details. The algorithm described here could also be used for the 
numerical evaluation of rational integrals. 

It was a pleasant surprise when, in the process of analyzing definite integrals 
of rational functions, we discovered that 

f°° cx i + dx 2 +e , 

t/fi = / — 7. 3 : — s dx 

J x e + ax 4 + bx 2 + 1 

admits a similar invariant transformation. We call this a rational Landen trans- 
formation. In the case of U§ the dynamical system (fTTj) is replaced with 

a n b n + ba n + bb n + 9 
~ (a„ + 6„ + 2)V3 ' ( 14 ) 
a„ + b n + 6 
" +1 " (a n + b n + 2)2/3' 

with similar rules for c„, d„, and e„. The derivation of (|14[) appears in [2]. 

The sequence (a n ,b n ) converges to (3,3) precisely for those initial data 
(ao,6o) for which the integral Uq is finite. Moreover, for the numerator pa- 
rameters, we have (c n , e?„, e n ) — * (1,2, 1)L, for some real L. The convergence 
of this method is discussed in [2], [12], and [14]. The invariance of Uq yields the 
identity 

* - k 

exactly as in (fT3"l) . Observe that ([6]) is also of this type: an integral given as 
the limit of an iterative process. Transformations similar to Q14[) have been 
produced in [3] for any even rational integrand. 

Until now all rational Landen transformations were restricted to even ratio- 
nal functions. In this paper we present the simplest example of a technique that 
we expect will extend to the general case (see [17] for details). 

The identities ([6]) and (|13[) yield iterative methods to evaluate the corre- 
sponding integrals. For example, the first four iterations of the evaluation of 



I = 

using © are given in Table 1. 



dx 



Ax 2 + 3x + 1 



4 



Table 1. 



n 


a n 


b n 


Cn, 





4 


3 


l 


1 


1.0731707317 


0.6585365853 


1.7317073171 


2 


1.3322738087 


0.0186646386 


1.31360991700 


3 


1.3228754233 


4.644065 xlO" 7 


1.3228758877 


4 


1.3228756555 


7.154295 xlO" 21 


1.3228756555 



The example presented in the table exhibits cubic convergence, faster than 
the convergence of the AGM. The exact value of / is 2n/V7, and ([6]) yields 
lim a n — \/7/2. The reader can check that the value 04 gives y/7/2 correct to 

n — >oo 

ten digits of accuracy. 

At the end of the amazing numerical calculation that led him to establish 
the invariance for the elliptic integral G(a,6), Gauss commented in his diary 
that this will surely open up a whole new field of analysis. This statement 
is certainly true. The reader will find in [5] a detailed discussion of how the 
arithmetic-geometric mean plays a fundamental role in modern computations 
of the digits of tt. This technique has also been used in [4] to create new and 
efficient methods to evaluate elementary functions. 

Over the years many proofs of (fT2| have been discovered. A number of them 
can be found in [18]. The authors are particularly fond of the succinct proof by 
D. J. Newman [19]: use x = b tan 9 and follow with x y~ > x + \J x 2 + ah. Change 
of variables is an art. 



3. THE QUADRATIC CASE. The goal of this section is to present the 
algebraic techniques that produce the transformation ([3]). We scale the inte- 
grand by multiplying both the numerator and denominator by an appropriate 
polynomial. This scaling is one of the main ingredients in the formulation of the 
Landen transformations. The other one will be discussed in the next section. 
We are motivated by the identities 

sin(30) T ,/ „s cos(3#) 
cos cos J 

where 

U(x) = x 3 - 3x , V(x) = 3x 2 - 1. 
The task is to find coefficients Zo, z±, Z2, 23, Z4 and eo, ei, ei such that 

{ax 2 + bx + c)(zox 4 + zix 3 + Z2X 2 + z$x + Z4) (15) 
can be written as 

e U 2 {x)+ ei U{x)V{x)+e 2 V 2 {x) (16) 



5 



with unknown coefficients z% and ei that are functions of the original parameters 
a, 6, and c. There is so much freedom, it can't be hard. 

Matching (fTS)) with (fTl)]) yields a system of seven equations for the eight 
unknowns. We use the first five to solve for the coefficients Zj in terms of a, 6, c, 
and the e,-. To start, comparison of the constant term in (|15p and (|16p gives 

z 4 = c _1 e 2 . (17) 

Using this value, we find that the first-order coefficient Z3 satisfies CZ3 — 3ei + 
bc~ 1 e2 — 0, which yields 

z 3 = c- 2 (3ce 1 -be 2 ). (18) 

The next powers produce 

Z2 — c~ 3 (9c 2 eo — 3bce\ + b 2 e2 — ace2 — 6c 2 e2) 

and 

zi = c~ (—96c eo + 36 cei — 3ac e\ — 10c e\ — b e 2 + 2abce2 + 66c 62), 

respectively. Finally, 

z = c" 5 (96 2 c 2 e - 9ac 3 e - 6c 4 e - 36 3 cei + 6a6c 2 ei + 106c 3 ei + 6 4 e 2 
— 3a6 2 ce2 + a 2 c 2 e2 — 66 2 c 2 e2 + 6ac 3 e2 + 9c 62). 

This leaves the two equations that arise from the two highest powers, which 
we use to find the parameters e^ . We solve the x 5 equation for e2 in terms of the 
parameters a, 6, c, ei, and eo- Substituting this information into the equation 
for the leading term produces 

6(6 2 -3(a-c) 2 )e = a(36 2 - (a + 3c) 2 )ei. (19) 

The system has one degree of freedom, which we exploit to ensure that the 
Zi and ei are polynomials in the parameters a, 6, and c. We initially choose 
e = a{(a + 3c) 2 - 36 2 ), from which it follows that e x = -6(6 2 - 3(a - c) 2 ). This 
in turn yields e2 = — c(36 2 — (3a + c) 2 ). 

The expressions for the coefficients z.; reduce to the following: 

z = (a + 3c) 2 - 36 2 

zi = 86(0 — 3c) 

z 2 = -6a 2 + 106 2 + 44ac - 6c 2 

Z3 = 86(c — 3a) 

z 4 = (3a + c) 2 - 36 2 

and, just to reiterate, 

e = a((a + 3c) 2 - 36 2 ) 

ei = 6(3(a - c) 2 - 6 2 ) 

e 2 = c((3a + c) 2 - 36 2 ). 
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In the latter formulas we already see a semblance of the iteration (|3J . 



4. ENTER TRIGONOMETRY. In this section we complete the construc- 
tion of the Landen transformation and establish the invariance of the integral |T]) 
under it. We establish the vanishing of a special class of integrals that appear 
as intermediate steps in this construction. 

We start with ([1]) and use the change of variables x — tan 9 to produce 

/ = r /2 — 2 d ± -. (20) 

J-tt/2 a sin" 9 + b sin 9 cos 9 + ccos 2 9 

The identities 

■< „ n sin(30) „ o „ cos(30) 

tan 3 6>- 3tan6> = 3tan 2 6»-l = 

cos 9 cos 6 9 

that were the reason behind the choices for U and V are then used to obtain 
r /2 sin k 9 cos 4 - fc 9 d9 , . 



fc=0 



- n / 2 e sin 2 (36») + e x sin(30) cos(36>) + e 2 cos 2 (36») 



from the integral {I) after it has been scaled according to the procedure de- 
scribed in section 3. 

The elementary identities 



cos 4 (9 = | cos(46>) + § cos(2(9) + § (22) 

cos 3 9 sin 9 = \ sfn(40) + \ sin(20) 

cos 2 sin 2 = | - | cos(40) 

cos sin 3 = i sin(26 l ) - \ sin(40) 

sin 4 6» = | cos(40) - | cos(20) + | 

transform the expression for I to a linear combination of 



Sh = r' 2 ^ k °) d ° (k=24) 

J-n/2 e sin 2 (30) + e x sin(30) cos(30) + e 2 cos 2 (30) 



and 



tt/2 



C fc = / , C ° Smde — (* = (), 2, 4) 

J-tt/2 eo sin (36*) + e\ sin(30) cos(36 l ) + e 2 cos 2 (36 l ) 

The magic of the Landen transformations comes from the vanishing of many 
of these integrals. This reduces (|2"Tj) to an integral of the type ([20]) with new 
coefficients, resulting in the transformation rule ([3]). Indeed, for even k the 
integrals Sk and Ck vanish if k is not a multiple of 3. To verify this, replace 9 with 
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u = 9 + 7r in the definition of Sk- Using sin(fc[« — 7r]) = (— l) k sin(fcu) = sin(fcu) 
and cos(fc[u — n]) = (— l) fc cos(ku) = cos(ku), we arrive at 



S k = 



3,7/2 sin(fcu) du 



tt/2 eo sin (3u) + e\ sin(3u) cos(3u) + e 2 cos 2 (3u) 



Adding this to the original Sk and taking advantage of the periodicity of the 
integrand we get 

If 271 sm{ku)du 

Sk — — 



2 Jo eo sin (3it) + e\ sin(3u) cos(3u) + e 2 cos 2 (3u) 

Now, we observe that both sin(3u) and cos(3u) are invariant under shifts by 
27r/3 and 47r/3, so 

f 27T sin(fcu) + sm(ku - 27rfc/3) + sin(A:u - 47rfc/3) 
Jo e sin 2 (3u) + ei sin(3w) cos(3u) + e 2 cos 2 (3u) 

The numerator in the integrand is the imaginary part of 

^iku _|_ ^i{ku — 27T&/3) _j_ gj(feti— 47rfe/3) ^iku /-^ ^ — 2irki/3 g— 47rfei/3^ 

and the last sum is 3 or depending on whether 3 divides fc or not. 

We conclude that the only terms that contribute to (f2"Tj) are the constants 
in (HU). Therefore 

1 f 27r 3z 4 + ^2 + 3z 

1 = / y dU, 

16 j eo sin (3u) + e\ sin(3u) cos(3u) + e 2 cos 2 (3u) 

where we have again appealed to periodicity to extend the integral to [0, 2vr]. 
The change of variables 9 — 2>u leads to 

/ = i r 2 *** d9 , 

8 J-tt/2 eo sin 6 + e% sin 6> cos 9 + e 2 cos 2 

so we have returned to the original form (|20p but with different coefficients. The 
result in Q is obtained by using x = tan 9 and the following identities: 



3Z4 - 


f^2^ 


- 3z 




8ei 




3Z4 - 


fz2^ 


- 3z 




8e 2 




3Z4 - 


f -22 4 


- 3z 



(3a + c) 2 - 36 2 



= b 



(3a + 


c)(a + 


3c) - 


b 2 


3( 


a-c) 2 


-a 2 




(3a + 


c)(a + 


3c) - 


& 2 


(a 


+ 3c) 2 


-36 2 




(3a + 


c)(a + 


3c) - 


6 2 



(23) 



5. THE ANALYSIS OF CONVERGENCE. In the last two sections we 
have shown the invariance of |T]) under the Landen transformation ([3]). We now 
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conclude by establishing the convergence of its iterates as in ([5]). In particular, 
we show that the error 

e„ := (o„ - \\J Aac - b 2 , b n , c„ - ~-y/ 4ac - b 2 ) 

satisfies e„ — > as n — > oo. Moreover, we demonstrate cubic convergence: 

IK + i|| < C||e„|| 3 (24) 

for some positive constant C . 

The analysis of convergence is simpler in the variables x = a + c, y = b, and 
z = a — c. The dynamical system (|3]) translates to 



%n+l 

z n +\ 

Vn+1 



Ax n 3z 2 t 3y n 



^1™2 — ~.2 — ~2 



(25) 



Arp2 — ~i2 — ~2 

o 2 _ 2 

Aqr^ — 7/2 — ^2 

^•^n tin 



with initial conditions x$ = x, y$ = y, and zo = z. 
We now prove that 



n — >oc n — >oc 



lim x n = \/ x 1 — y 2 — z 2 , lim y n = lim z n — 

n — >oo 

or, equivalently, that 



lim ( x n - \J x 2 -y 2 - 



yl+zl = 0. 



(26) 



This is equivalent to (O, so it will finish the proof of convergence. 

To complete the change of variables we use the invariance of the discriminant 
(0 to obtain 

222 222,* i2 

and we write w = V Aac — b 2 . The first equation of iteration (|25p becomes 

' xi + 3w 2 ~ 



x n +i 



3x 2 + w 2 



(27) 



with initial conditions xq = a + c (> 0). The required limit in (|26|) is now 

lim x n (x n - w) = 0, (28) 

n — >oo 

with y n and z n absorbed into the constant w. The number of variables has been 
reduced from three to one. 

Using one last change of variables, q n — —ix n /w, we reduce (|27|) to 
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In ~ 3g» U{q n ) flyQ s 

What we need to prove in order to establish (|2"5j) is that q n — > —i, (The 
polynomials U and V introduced in section 3 have miraculously reappeared!) 
The trigonometric identity 

U(cot9) 

vwi) = cot(30) < 

coupled with a representation of the initial condition as 

q a = cot (it) (30) 
for some t (0 < t < 00), shows that (|29|) simplifies to 

E/(cot it) . , 

91 = vioT^ = cot(3ii) 

and, in general, 

e 2 * 3 " + 1 

q n = cot(3"it) = 



e 2t3" _ ! 

We conclude that q n — > — i, whence £ n — > w, as desired. 

To verify (|30p. we write go = — *^ with e? = (a + c)/y/4ac — b 2 . Now recall 
that 4ac — b 2 > 0. An elementary argument shows that d > 1, so we can take 

, 1 , „ 1 , d + 1 
t = coth x (d) = -ln-7— p. 

The fact that the convergence is cubic follows directly from 

2 



which decreases to like e -2 * 3 . This implies (|24[) and completes the proof of 
convergence. 
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